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LIMIT BEHAVIOUR OF THE TRUNCATED PATHWISE 
FOURIER-TRANSFORMATION OF LEVY-DRIVEN CARMA 
PROCESSES FOR NON-EQUIDISTANT DISCRETE TIME 

OBSERVATIONS 

ZYWILLA FECHNER AND ROBERT STELZER ^ 


Abstract. This paper considers a continuous time analogue of the classical autoregres¬ 
sive moving average processes, Levy-driven CARMA processes. First we describe limiting 
properties of the periodogram by means of the so-called truncated Fourier transform if 
observations are available continuously. The obtained results are in accordance with their 
counterparts from the discrete-time case. Then we discuss the numerical approximation 
of the truncated Fourier transform based on non-equidistant high frequency data. In 
order to ensure convergence of the numerical approximation to the true value of the 
truncated Fourier transform a certain control on the maximal distance between obser¬ 
vations and the number of observations is needed. We obtain both convergence to the 
continuous time quantity and asymptotic normality under a high-frequency infinite time 
horizon limit. 


1. Introduction 

The classical autoregressive moving average process ARMA has been broadly discussed 
in the literature. For a comprehensive discussion see e.g. the monograph by Brockwell 
and Davis and references therein. In the discrete time models we restrict ourselves to 
observations at hxed equidistant points in time. In many cases these observations made 
at discrete times come from an underlying continuous process, thus the natural question 
arises: can we model also the time series in continuous time? One of the earliest results 
dealing with properties of such processes can be found in Doob |11] . Later this problem 
was discussed by Brockwell in for continuous time ARMA processes driven by Gaussian 
noise. The next step was to extend these ideas to the models with noise modelled by 
jump processes, so-called Levy-driven CARMA models introduced by Brockwell in j^. In 
these papers time series are modelled as continuous time processes with continuous time 
noises (with or without jumps) and the inference is based mainly on discrete equidistant 
data. One of the latest results can be found in the paper of Brockwell, Davis and Yang, 
which consideres QML estimations of the AR and MA parameters based on equidistant 
observations. 

The estimation procedure of Levy-driven CARMA processes in high-frequency settings has 
been discussed by Fasen and Fuchs in |13] , where the authors deal with the limit behaviour 
of the periodogram of CARMA processes under equidistant sampling when the sampling 
interval tends to 0. The results are analogous to ARMA processes: the periodogram 
for CARMA processes is not a consistent estimator of the spectral density function, but 
after appropriate smoothing the consistency can be obtained. Some related results were 
discussed by Fasen and Fuchs in (14) , where asymptotic distributions of periodograms of 
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CARMA processes driven by a symmetric a-stable Levy noise are obtained and where it 
is shown that the vector composed of periodograms for various frequencies converges in 
distribution to a function of a multidimensional stable random vector. Likewise, Fasen 
considers the behaviour of the periodogram for an equidistantly sampled continuous time 
moving average process when only the number of observations goes to inhnity. 


The problem of statistical analysis of such processes has been studied further for example 
by Gillberg in his dissertation (I^ , where different approaches to the estimation of CARMA 
processes with Gaussian noise are discussed both using equidistant and non-equidistant ob¬ 
servations. The author works mainly in the frequency domain. He describes the properties 
of the truncated Fourier transform of a CARMA process with Gaussian noise on a hxed 
interval [0, T] based on equidistant frequencies. In the non-equidistant case he has used a 
method based on splines in order to hnd an approximation of the spectral density. 


Another approach for the estimation of a zero-mean stationary process (Yt)t£M. with hnite 
second-order moments and continuous covariance function has been discussed by Lii and 
and 17 , where they described some properties of a smoothed periodogram. 


Masry in 

Here observations are assumed to be given on a random grid (rfc) of an interval [0,T], 
where is a stationary point process on the real line which is independent of 


In the present paper we are going to describe the asymptotic behaviour of the so-called 
truncated Fourier transform of a CARMA process, which is a building block for an estima¬ 
tion of the spectral density of a CARMA process. We are going to use some of the ideas 
from |15| to prove results in more general settings. 


The paper is structured as follows: hrst we recall second order Levy-driven CARMA mod¬ 
els and summarize the results needed later in Section |2] Then we dehne in Section |3]the 
truncated Fourier transform of a CARMA process and we investigate its asymptotic prop¬ 
erties at a hxed frequency: for a non-zero frequency we obtain that the limiting law of the 
real and imaginary part is the two dimensional normal distribution with mean zero and the 
covariance matrix depending on the spectral density of the CARMA process. If we con¬ 
sider the truncated Fourier transform at zero, we obtain a one dimensional normal law with 
mean zero and variance depending only on two parameters of the CARMA process. We 
show that the limiting law of the joint distribution of the squared modulus of the truncated 
Fourier transform at different positive frequencies converges to a vector of independent and 
exponentially distributed random variables with mean depending on the values of the spec¬ 
tral density. All these results can be interpreted as the limiting behaviour of the truncated 
Fourier transform when the CARMA process is observed continuously. The next step in 
Section 3.2 is to approximate the truncated Fourier transform when the CARMA process 
is observed on a non-equidistant deterministic grid. In order to hnd a numerical approx¬ 
imation value of the truncated Fourier transform we apply the trapezoidal rule. We are 
interested in the convergence of the truncated Fourier transform when the length of the 
interval T goes to inhnity and the mesh of the grid to zero. Since the interplay of the 
length of the interval, of the number of elements of the grid and of the maximal distance 
between the elements of the grid plays a crucial role, in order to ensure the convergence 
of the approximating sum to the true value of the truncated Fourier transform we have 
to impose some limiting conditions on these quantities. In the last Section]^ we look at 
some illustrative simulations of the truncated Fourier transform based on non-equidistant 
observations. We consider Ornstein-Uhlenbeck type (CAR(l)) and CARMA(2,1) processes 
driven by a standard Brownian motion, a Variance Gamma process and a “two-sided Pois¬ 
son process” and we compare our simulations with the theoretical asymptotic distributions 
described earlier. 
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Notation. The symbol N := {1,2,3,...} denotes the set of positive integers. No := 
N U (Oj, M is the set of real nnmbers and C denotes the set of complex nnmbers. The 
symbol resp. denotes the space of real- (resp. complex-) valned matrices 

with m rows and n colnmns. For A G the symbol denotes the transposed of a 

matrix A. We are working on a given filtered probability space (11, T, IP) satisfying 

the nsnal hypothesis (cf. Protter |20] , Chapter 1). 

Moreover, hy X = Y we denote that the random variables X and Y are eqnal in distribn- 
tion. 


2. PRELIMINARIES 

We begin with the model set-np given by Brockwell (see [^, |^). A second-order Levy- 
driven continnons-time ARMA(p, q) process is defined in terms of a state-space represen¬ 
tation of the formal differential eqnation 

(1) a{D)Y{t) = b{D)DL{t), t > 0. 

Here, D denotes differentiation with respect to t, non-negative integers p, q satisfying 
p > q and is a one dimensional Levy process (i.e. a continnons time process 

with stationary and independent increments and L(0) = 0 a.s.) with EL(1)^ < oo. A 
comprehensive monograph dealing with Levy processes is e.g. [^. The polynomials 

a{z) := -h aizP~^ H-h Op, b{z) := bo + biz-\ -h bp_iz^~^, 

are called the autoregressive- and moving average polynomial, respectively. We assnme 
that bq ^ 0 and bj = 0 for q < j < P- The state-space representation consists of the 
observation and state equations: 

(2) Y{t) = b^X(t), 

(3) dX{t) = AX{t)dt + edL{t), 

where 



0 

1 

0 

0 ■ 


■ x{t) - 


0 

0 

1 

0 


XW(t) 

A := 

0 

0 

0 

1 

, X(t):= 

A(P-2)(t) 



ttp—1 

ttp—2 • • ■ 

-ai_ 


x(p-i)(t)_ 


e := [0,... 0,1]^ , b := [6o, 6i,... bp-if , 


i.e. 

X{t)eMP^^ eeRP, heMP. 

If p = 1, we set A = —ai. 

Assumption 2.1. The Levy process satisfies EL(1) = 0 and E|L(l)p = < oo. 

Observe that K[L{s)L{t)] = min{s, t}E|L(l)p. It was shown by Brockwell in that the 
solution X(t) of (|^ satisfies 

(4) X(t) = e^*X(0) + [ e^(‘-“)e(iL(n), 

Jo 

where the integral is defined as the L^-limit of approximating Riemann snms. 
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Assumption 2.2. X(0) is independent of {L{t))t>o. 


From now on let us assume that Assumption |2.2|holds. It is well-known ([^ Proposition 2]) 


that under Assumptions 2.1 and|2.2|the process {X(t)}i>o is strictly stationary and causal 


iff X(0) has the same distribution as e^'^edL{u) and the p (not necessarily distinct) 
eigenvalues Ai,..., Ap of A have negative real parts, i.e. 

5R(Ai) <0, i = l,...,p. 

Now we extend the Levy process {L{u))u>o to the whole line in the usual way: Let L = 
{L{t))t>o be an independent copy of {L{t))t>o- For t G M we define 

L*{t) := L(t)l[o^oo)(t) + -L(—t—)l(_oo,o](t)- 
In order to get stationary solutions of ^ we need the following assumptions: 
Assumption 2.3. All eigenvalues of A have strictly negative real parts. 

Assumption 2.4. 


X(0) = f 
J — c 


e ^'^edL*{u) 


In Brockwell it was shown that if Assumptions 2.3 and 2.4 are satisfied the process 
{X(t)}tgR given by 

(5) X(t) = f 

J —OO 

is a strictly stationary solution of ^ (with L replaced by L*) for t G M with corresponding 
CARMA process 

(6) Y{t) = j b'^e^(*-“)edL*(n). 

J —OO 

For t > 0 one can rewrite it in the following form 

(7) y(t) = b^e^‘X(0) + / b^e^(‘-“)edL(n). 

Jo 

In the present paper the spectral density of a CARMA process plays a crucial role. 
The spectral density is the Fourier transform of the autocovariance function 'Jyih) : = 
Cov{Y{0),Y{h)) for /i G M. The spectral density of a CARMA process is 


( 8 ) frita) = 


1 


-fY{h)e-^^dh = 


cr^ \b{iuj)\" 
2 tt |a(*Ci;)| 


2 ’ 


w G 


3. Limit behaviour of the Fourier transform 

In this section we are going to deal with the Fourier transform of the CARMA process 
assuming that the observations are given continuously on the time interval [0, T]. A similar 
idea for Gaussian CARMA processes was presented in |15] for equidistant observations. 
The truncated continuous-time Fourier transform of the process T at a fixed frequency 
a; G M is given by 

T-tRXu,) := Y 

Observe that the norming constant is taken as this ensures convergence in distribution 
for T —7- OO as will be shown later. 
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3.1. Properties of the Truncated Fourier Transform of a CARMA Process. First 
we derive an alternative representation. 

Lemma 3.1. Let X and Y be proeesses given by the state-spaee representation (|^ and 
(|^. Suppose that Assumptions 2.1 2.2 and 2.3 are satisfied. Then the truncated Fourier 
transform of the CARMA process Y at a fixed frequency uj is of the form 

(9) 7-t(F)(l,) = ^ ~ ' 

or equivalently 

(10) YT{Y){ov)~h^{iuI-A)-^ 

edL{u) + (l - X(0) 


Proof. Let UJ be an arbitrary frequency. Observe that by Corollary 3.4 from 21 p. 51] 
one has 

b{iuj) 


(A — ioji) ^e =-— 


a{ioj) 


Denote 


F{t) = b^(A - t e [0, T], 


G{t)= / e-^'^edL{u) t e [0,T]. 

Jo 

Observe that G(0) = 0 and since F is continuous and of finite variation, we get [F, G] = 
0, where [•,•] denotes the usual quadratic covariation of semimartingales (see e.g. (^). 
Applying the (multidimensional) integration by parts formula 

[ dF{t)G{t) = F{T)G{T) - F(0)G(0) - [ F{t)dG{t) - [F, G] 

Jo Jo 

= F{T)G{T)- [ F{t)dG{t) 


we obtain 


[ dF{t)G{t)= [ h'^{A - iujI)-\A - iojI)e^^-^^^J [ e-^^edL{u)dt 
Jo Jo Jo 

= [ f h^e^J-'^'^edL{u)e-^^^dt 
Jo Jo 

= b'^(A — iujl)~^e 


TtK .r\-l„(A-ia;/)T f „-At f .T(x_- T--.-l (A-iujI)t -At 


e-^^edL{t) - b^ (A - iujI)-^e^^-^‘^^>^e-^^edL{t) 
10 Jo 


= W {A — iuji) ^e 




eMT-t)^dL{t) + 


r ^-iLot 
a{ioj) Jq 


dL{t). 


Thus 

( 11 ) 


fo ^^edL(u)e ^^^dt 

= b'^(A - e^(^-*)edL(t) + ^ e-^‘^*dL(t). 
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Using the form of the strictly stationary solntion of (|^ given in Q we get 

(12) [ e^^'^-^^edL{t) = X{T) - e^^X(O). 

Jo 

Moreover, since = {iuil — A)~^{I — have 

(13) [ h^e^^-^‘^^Jx{0)dt = h'^iiujl - A)-^{I - e(^“*‘^-^)^)X(0). 

Jo 

We have 

^ A r (b’'e^‘X(0) + j‘ b’'e^(‘-“>eiiL(u)) e-‘“‘dt 

0.0.0 + yf- e'‘“’'x(r)). 

To get the eqnivalent form note, 


VfF{Y){u) = h^{iujl - A)-^ 


e-^^'^dL{u) + (X(0) - e-*‘^^X(r)) 


112 


^ [iujl - A) 


-1 


^ _ g-ib.TgA(r-«)^ edL{u) +{l- X(0) 


which completes the proof of this Lemma. □ 

The next step is to calcnlate moments of the trnncated Fonrier transform. First, recall the 
so-called compensation formula: If {Lt)t>o is a Levy process with hnite hrst moments and 
/ is a bonnded deterministic fnnction, then 


(14) E 


L4o 


f{u)dL^ 


= E[Li] / f{s)ds. 
Jo 


Secondly, observe that the solntion of the system ([^ and (|^ is of the form where X 
is the process with mean m{t) = E[X(t)] and Px{t) = iE[X(t)X(t)'^] satisfying 

mx{t) = e^*mx(0) 

(15) Px{t) = e^^Px{0)e^^^ + f du 

Jo 

In particnlar, for stationary processes these solutions are constant and the so called Lya- 
pnnov eqnation 

(16) APx + PxA"’" + cr^ee^ = 0 


holds true. For Levy-driven CARMA processes the form of the antocovariance fnnction in 
terms of solutions of Lyapnnov equations is formnlated e.g. in 19, Proposition 3.13.]. 


We are hrst going to show that the trnncated Fonrier transform of a stationary CARMA 
process is a zero-mean random variable. Next, we hnd the covariance between the trnncated 
Fonrier transform at two different freqnencies. As we have mentioned earlier, the spectral 
density fnnction plays a central role. 
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Theorem 3.2. Let X and Y be processes given by the state-space representation and 
(|^. Suppose that Assumptions 2.1. 2.2 2.3 and 2.4 are satisfied. Then E(J^ 7 ’(y)(a;)) = 0 
for all a; G M. For uJi,U 2 G M we have 

(17) E[J^t{Y){u;i)J^t{Y){uj2)]=ct^Pp^^ + ^K{T,u;i,-ui), if =-U2 


|a(*wi)|^ 


and 


(18) E[TTiY){coi)J^T{Y){uj2)] =-KiiT,oji,co2), if coi ^-C 02 , 
where K is a bounded function of T given by ( |20[ ) below and 

Ki{T,uJi,ui2) = K{T,LOi,LV2)+h'^{ioJiI-A)~^a‘^- — ^^(^1 + ee'^(zcj27-A'^)~^b. 

^(a;l + UJ 2 ) 

Proof. For the first part it is enough to observe that by the compensation formula 


E 


T 


Q e ^‘^'^dL{u)j = 0 and E[X(t)] = 0. For the second part observe that using Lemma 
3.1 and formula (10) we have 

E[TT{Y){iOi)TT{Y){u2)] = ^h^{iiOiI-A)-^x 


E 


_ ^-iuiiT^A(T-u) 


qT f^-io^2U _ g-jaJ2TgA^(T-M) 


1 


) edL(u) + (^ - e(-^‘^Y+A)T^ X(0))^ 
dL(u) + X(O)'^ - e(-^‘^Y+Ar)T^ 


{ioj 2 l — = —b’^(it<;il — A) ^I{iuj 2 l — AJ') ^b, 


where I = Ii + I 2 + H Y Ia with 

pT 

Ii := E ' 


I 2 :=E 
h :=E 

/4 := E 


/o 


/o 


gT t^-iuJ2U _ ^-iuJ2T(T-u) 


^AiT-u)^ edL{u) 

_ g-io.irgA(r-«)^ edL{u) ■ X(0)^ {l - eY^ 2 i+Ar')T\^ 

/.T 

I _ g(-*‘^i^+A)r') 

J _ Q{-i‘^\I+A)T\ 


J X(0) • 


'0 


gT (^-iuj2U _ ^-iuJ2T^A'T(T-u) 


dL{u) 


'0 


dL{u) 


) 


X(0)X(0)^ [I-e 


^-1^.021+AT )T 


We have that I 2 = I 3 = 0 since {Lt)t>o is independent of X(0). Observe that by the ltd 
isometry, the compensation formula and the fact that E[[L, L]i] = Var(L(l)) = we have 


It :=E 


e-*‘^i“edL(n) / e'^e-*‘^2“dL(n) 


= E 


e-^^‘^^^‘^^^^ee^d[L,L], 


pT pT 

= E[[L,L]i] / = du. 

Jo Jo 


Thus 


J a^Tee^, =-^ 2 , 

(iy) — A 2 l~6xp(—Ti(aJi+t.t;2)) T / 

r - iL+L) a;i^-a;2. 


i(aJi+W2) 

Thus, if oji = — 0 J 2 , then 
1 


—b^ (iwil — A) ^l{{ioj 2 l — A^) = — ■ a^Th^ {iujil — A) {ioj 2 l — A 


,r\-i 


b 
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= (A — iojil) {iuiil + 

_ J.2 f ( b{-iuji) \ _ 2 l^(^^i)l 


a{iuji) J \ a{—iuji) J |a(iu;i)P 


Now 


It :=E 


Jo 


= e“*‘^2T]£ 




= e-^‘^^^E[[L, L]i] / 
Jo 

Jo 

In the same way 


u 


It ■.=E 


Ljo 


e-*‘^i^eA(^-“)eciL(u) / e^e-^^^^dL{u) 


Jo 

Combining these two we arrive at 

if + if =e-*(‘^i+‘^2)r^2 ee'^(ia;i/ + A^)"^ J'g(*a;i7+A^)T _ 


+ {iuj2l + A)-^ L{io^2l+A)T _ J 


Now 


if ; = E 


ee 


^ f ^T-iuj2TA'^{T-u\ 


Llo 


>edL{u) / 
Jo 


'>dL{ 


u 




Now 


I4 =E[X(0)X(0)'^] - e"*‘^i^e^'^E[X(0)X(0)^] - e“*‘^2T]gjx(o)X(0)'^] 

+ e“*(‘^i+‘^2)rgATj£^X('Q^X('Q^rjgArr^ 

By stationarity we have 

E[X(0)X(0f ] =: Px = Px{0) = Px{T), 
where Px satisfies (16). Combining this with (15) we obtain 

if + I 4 =Px - e-^^^^e^^Px - e“*‘^ 2 Tp^gA^r ^ ^at 

+ e-*(‘^i+‘^2)r(PY _ e^^Pxe^^^) 


^g-iojiTp^ (I - + e“*‘^2T 


I-e^^)Px 


+ Px{l- - e- + e 


^-iuJiT _ ^-iuJ2T _|_ ^-i{uJi+uj2)T 

Since A is a stable matrix, is bonnded 
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Thus 


( 20 ) 


K(T,lji,lj 2) =h^ {iujil — A) 


-1 




ee 


+ A^)-^ f^^{imi+AT)T _ j 


+ {iU 2 l + A)-^ L{i^ 2 l+A)T _ J 


ee 


+ e” 


I-e 


+ e 


—iu!2T 




+ Px(l- g-i(a;i+a; 2 )r^ J 

is bounded in T for fixed a;i,a ;2 G M. □ 

Now we give the form of the covariance matrix. Put 


S(a;i, W 2 ) 


[5jij]l<ij<4 — IE 


'^Pt{Y){o)i) 


'iflPT{Y){u)i) 

r-| 

^Pt{Y){o)i) 


^Pt{Y) (ooi) 


^Pt{Y){o)2) 


iliPT{Y){0)2) 


<APt{Y){o)2) 


^Pt{Y){0)2) 



Theorem 3.3. Let X and Y be processes given by the state-space representation ([^ and 


(|^. Suppose that Assumptions 2.1 2.2 2.3 and 2.4 are satisfied. For oji 002 and 


LOi — 0)2 there exists a bounded matrix K 2 G such that 


T.{u)i,U)2) 


^ 2.( \h{iU)l)? \b{iu),)\^ \b{iu)2)? \b{zu)2r \ 1 

2 ^ |a(iw 2 )| 2 ’ \a{iu2W) T ^ 


Proof. For k,l = 1,2 let us denote 

Si(wi,a; 2 ) ■.= E[iRPT{Y){u)i)iRPT{Y){u) 2 )], ^2(0)1,0)2) := E[<APT{Y)iu)i)<APT{Y){u) 2 )], 

Es{u)i,u)2) := E[3f?TT(I^)(wi)9.FT(i^)(a;2)] • 

All entries Sjj of the matrix S are of one of the above forms. Indeed, Sn, S 33 are of the 
form Si for k = I and k, I G {1, 2}. Similarly, S 22 , S 44 are of the form S 2 for k = I and 
k,l G {1,2}. Moreover, S 13 , S 31 are of the form Si for k 7 ^ I and k,l G {1,2} and S 24 , 
S 42 are of the form S 2 for k I and k,l G {1,2}. All other elements are of the form S 3 . 


Observe that for each u) we have 

Using Theorem |3.2|we obtain 


QPt(Y){u)) = 


Pt{Y){u)) - Pt{Y){-u)) 


2i 


Ei{u)l,U)2) 


: = 


S 2 (a;i, ^ 2 ) 




2 ^ 


2 |fe{ 0 ) 

2 

|a( 0 ) 

T 

b{uji) 

2 


2 


2 

aUiJi) 

2 


^ + T^( 0 ), 

2 I 

^ + 7pAi^i(wi), 

2 . 

2 + 7 pAi^ 2 (<Vl), 

0 , 

k'ISuF " T^ 2 , 2 (wi), 


|a(jtJi)P 

y A 2 , 3 ( 1 ^ 1 , a; 2 ). 


0Ji= u)2 = 0 ; 

U)l = OJ 2 ] 

CJi = -u)2-, 

u)i 7 ^ UI 2 andwi 7 ^ — 0 J 2 , 

u)i = 0 or u )2 = 0 ; 

UJl = 0)2] 

CJI = -U)2; 

u)i 0)2 and cji 7 ^ —(^ 2 , 


S3(a;i,a;2) 


0 , 

^A3^i(wi), 

7 A 3 , 2 (^ 1 , W 2 ), 


0)2 = 0 ; 

0)1 = 0)2 or 0)1 = — 0 ) 2 ; 
cji 7 ^ 0)2 and cji 7 ^ —a; 2 . 
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Here K is given by (20) and Kij are bonnded in T for i,j = 1,2,3. □ 

Now we are going to investigate asymptotic properties of the trnncated Fonrier transform. 
First, we will show that the second snmmand of ([^ converges in probability to zero. 

Lemma 3.4. Let X and Y be proeesses given by the state-spaee representation ([^ and 
([^. Suppose that Assumptions 2.1\ 2.2 and 2.3 are satisfied. Let 


Z{T) ■.= TTiY)iu)- 


1 b{iLo) 

Jq 


f 


—iut 


dL{t). 


Then 


- hm \Z{T)\=t). 

T^oo 


Proof. Observe that 


|Z(T)| = 


1 


[X(0) - e--^X(r)] 


Obvionsly, 


< ^ \h^{iu;I - A)-^X{0)\ + ^ \h^iicvL - A)-^XiT)\. 
lim —^ \h^(icjl — j4)“^X(0)| = 0 a.s. as T —)> oo. 

T->-oo 


Becanse of stationarity, X(T) is bonnded in probability and converges to zero thns 


-^\h^{iuL-A)-^MT)\ 


0 in probability. 


Therefore 


- hm |Z(r)|= 0 . 
T^oo 


This completes the proof. □ 

Now we will show that the first snmmand of formnla (|^ converges in distribution. Thns, 
together with Lemma 3.4 we obtain the limit in distribntion of the trnncated Fonrier 


transform. We have two cases: the first case is if the freqnency w = 0. Then the trnncated 
Fonrier transform is a real valned function. In the second case for freqnencies a; 7 ^ 0 
the trnncated Fonrier transform is a complex valned fnnction. In both cases we first give 
the description of the distribntion of the trnncated Fonrier transform and afterwards we 
describe the distribntion of the sqnared modnlns of the trnncated Fonrier transform. 

Theorem 3.5. Let X and Y be processes given by the state-spaee representation ([^ and 
([^. Suppose that Assumptions\2. j\ and 2.3 are satisfied. Let 


Z{T) : = 


J_M f 

\/ra(0) Jq 


dL{t). 


Then 


d - lim Z(T) ~ 7VA 0, 

T^oo \ 


m 

a(0) 


a 


d — lim —^ 
T^oo 


a{Q)Z{T) 


6 ( 0 ) 


X^(l) 


Proof. Observe that dL(t) = L{T), thns Z{T) = By the standard Central 

Limit Theorem d — lim'r_,.oo -^L{T) = Af{0,a‘^). Therefore d — lim 7 ’_,.oo = 

^rfn 
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Observe that for all re G N the random variable ~ Af (0,1). Then by the continuous 

mapping theorem we have d — limr_>.oo ^ 

In order to find the asymptotic distribution of the truncated Fourier transform we use the 
multivariate Central Limit Theorem. Note that we state all results for positive frequencies 
as the corresponding results for negative can be obtained by taking the complex conjugate. 


a(0)Z(T) 

m 




Theorem 3.6. Let X and Y be processes given by the state-space representation ([^ and 
(|^. Suppose that Assumptions \2.i\ and \2.S\ are satisfied. Assume that w > 0. Put 

1 b{iuj) r 


Z{T) : = 


y/T a{iu) 


dL{t) 


and 


Then 


where S 


b{iLo) 

a(iuj) 


2 

hx2- 


Z{T) 


'^Z{T) 

QZ{T) 


d- lim Z{T) Af{0,Y), 

T—>-oo 


Proof. We firstshow that /o ® '^^dL{t) is asymptotically normal. For iV G N and 
j G {0,...,iV- 1} put 

- r2-K{j + l)/u] ^ 




R1 





f27Tjl(jJ 

2'K{j + l)/0J . 


cos{ujt)dL{t) 


r 

_J27rjlcj 


sm(ujt)dL(t) 


Observe that Xj are independent and identically distributed random vectors with mean 
zero and the covariance matrix Si := ^^/ 2 x 2 - Therefore, 


—iut 


Applying the classical CLT we obtain 


j=0 

!— 2 wN 

/o “ e-^^^dL{t) 


N-l 


So 


/2ttN jo 


-f 

Vn Jo 


N-l 

dL{t) = Y, Xj- 

j=0 


e-^‘^^dL{t) ^ W ~ W(0, Si) as N ^ oo. 


W(0,Si), where Si = ^Si = ^/ 2 x 2 - Put 


A := 





b(iuj) 

a(iij) 

b(iLj) 

a(iLj) 


Observe that 

A- 


yd c2-kN/lj 
'2-kN -iO 
yd I'2itN/lU 

Jc 


IVmx JO 


JZTT 1 \/U 

sm{ujt)dL{t) 



'^Z{2tiN) 


<AZ{2'kN) 


Thus Z = A - X IS normally distributed with mean zero and the covariance matrix S = 
2 


ASiA^ = ^ 


b{iu}) 

a{iu) 


h 


x2- 


□ 
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Now we apply this theorem to find the asymptotic distribution of the truncated Fourier 
transform squared. 

Theorem 3.7. Let X and Y be processes given by the state-space representation ([^ and 
IL Suppose that Assumptions\2. i\ are satisfied. Let Z be defined as in Theorem 


3.6 Then |Zp ~ Exp ( 


with mean A. 




a{iu)) 


Proof We use the notation of the proof of Theorem 


, where Exp(A) denotes the exponential distribution 

Thus |Zp is proportional to 


3.6 


chi-square random variables with two degrees of freedom, i.e. |Zp = ^ 


X~x 2 ( 2 ). Thus |Z|2~r 1,^ 




so |.Z’p ~ Exp ( fj 


b(iuj) 


b{iu}) 

a{iu) 

□ 


X, where 


Now we are going to give the description of the convergence of the random vector consisting 
of the truncated Fourier transform at different frequencies. 

Theorem 3.8. Let X and Y be processes given by the state-space representation ([^ and 
([^ . Suppose that Assumptions \2.1\ l2.2l and |.2.3| are satisfied. Let 0 < uji <■■■< be 

fixed frequencies. Then [5R ( 7 ' 7 ’(T)(a;j)), O'^ converges to A/" ^0, 
with 


B = diag 

iT 


b{iuji) 


a (foil) 


b{iuji) 


a{iuJi) 


b{iujd) 


a{iujd) 


b{iujd) 


a{iuid) 


and ^ converges to a random vector whose coordinates are independent 

Exp ( 




b{i<jjj) 


a(iLJj) 


distributed random variables for j = 1,... ,d. 


Proof For fixed n G N and k = 1,... ,n, put 

p2k7v p2k'K 


X 


[oji) := / cos{oJit)dL{t), (coi) := / sm(uJit)dL(t), i = l,...,d. 

J2(k—l)7r J2(k—l)7r 

2 . , . . . / /o.a \2 


Let (^Sn' Var i)(a;i)] and Var [x(2*)(a;j)]. 


Put 




ELiVf’G.) 


A 2 *-l) 


J2i) 


Then we will show that by the Cramer-Wold-device the random vector Z G 


T>2d 


with 


Z = 


X, 


(2i-l) 


Wj), 


1 T 


converges in distribution to JY ( 0 , 12dx2d)- 

We first apply the Lindeberg-Feller Central Limit Theorem (see e.g. Billingsley [^) to each 
coordinate of the vector Z. Observe that for all i = 1, ..., d by the Ito isometry we obtain 


Var (Xf = Var 


r*2/c7r 


r2k7r 


f 2{k—l)'K 


cos{uJit)dL{t) = 


' 2(fc— 1)77 


cos^ {ujit)dt 


= a 


2 4:TTCOi + sin(47rt<;iA:) — sin(47ra;i(A; — 1)) 
4:UJi 


A2i-l)l _ 


Var 


X 


(2i-l) 


(uJi) 


= a 


k=l 


2 4n7ra;i -|- sin(47ra;in) 
Auji 


Thus 
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In the same way, 






= a 


Observe that 


k=l 


lim - = lim - = 

1^00 n V / n^oo n \ J 


2^mTUJi — sin(47ra;in) 
4(^2 


1 / ^9,^2 


cr^TT. 


If the Lindeberg condition is satisfied, the 2z-th, respectively 2i — 1-th coordinate of Z for 
i = 1 ,..., d, i.e. 

f cos{uJit)dL{t) 


a^/AnnuJi + sin(47rnwi) Jo 

,r 

crydvrnwj — sin(47rnwij Jo 


sin{Ldit)dL{t) 


converges to AA(0,1). Taking 

limuji 

and noting that 


a^/AimoJi — s\Ti{ATmuJi) 


2\J2'KnijJi 


lim = lim ^ 

n^oo y2 n-5>oo y2 

is constant at all frequencies, by Slutsky arguments for i = 1, ..., d we get 

1 fSTTn / 2 \ 

(21) cos{uit)dL{t) = M ^0, yj , 

1 /■27rn / 2 \ 

(22) — sm{u:,t)dL{t) = z(^^\co,)yJ^^\co,) ^ Af (^0, -j . 

Now we are going to prove the Lindeberg condition for odd coordinates of Z (for the even 
ones an analogous reasoning holds), i.e. for all e > 0 it holds 


lim 


1 




n^oo / (21-1) , 

Sn k=l 






= 0 . 


Observe that if the random variables {xj^^ ^^(wj)} are uniformly square integrable, then 
they satisfy the Lindeberg condition. Indeed, 


42 i-l) 


-2 


E*: 


fc=i 






= ( 


-2 


E® 


fc=i 




< Sn 


_2 

( 2 *-l)'\ ^ g^p ]£ 


xt a 


-1 


2 dnvrcjj -I- sinUiruJin] . ^ 

= I a --- n sup E 

4Wj / fc=l,...,n 




{2i-l) 
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Since lim„^oo 4 n 7 rcj,+sm( 47 ra;,n) ^ ^ uniform square integrability implies in this 

case the Lindeberg condition. It remains to show the uniform square integrability of 

I ^ JfceN 

Assume first, that our driving process {L{t))t>o is of bounded variation. Then 


Mk = 


r2kn 


'2{fc —l)7r 


cos{LOit)dLi 


< 


p2k'K p2k'K 

/ \cos{uJit)\d\Lt\ < / d\Lt\, 

J 2{k — l)7r J 2(k—l)7r 


2k'K 
2{k—l)'K 


where | • | denotes the total variation of the process. But J 2 (k-i)n ^ We 

have 


E [\Mk\‘^l^\Mk\>K}] < E 
= E 



2 



2 

^{|/oW|L,||>R} 


By the square integrability of d\Lt\, which is implied by the square integrability of 
{L{t))t>o we obtain the uniform integrability of {Mk)k£n- 

Now we assume that {L{t)) is a square integrable martingale with finite moments of all 
orders. Observe that Xk is square integrable for all fe G N. By the Burkholder-Davis- 
Gundy Inequality (see e.g. Protter |20] ) for each p > 1 there exists a positive constant Cp 
such that 



/ i'2kn \ ^ 



p2kn 

p2k7T 

p/2- 

E 

I / COs{LOit)dL{t) J 

< CpE 


/ cos{oJit)dL{t), 

/ cos{uJit)dL{t) 



\J2{k-l)n J 



^2{fc —l)7r 

/2{k—l)7T 



Since 


p2k7T p2k7T 

/ cos{u}it)dL{t), / cos{ujit)dL{t) 

l2(k—l)7r ^2(fc —l)7r 


r*2/c7r 


/2(fc-l)7r 


cos^{uJit)d[L, L]t 


using the above inequality for p = 4 we obtain 



/ p2k7v \ 


p2kn 

E 

I / cos{uJit)dL{t) J 

< C 4 E 

/ cos^ {uJit)d[L,L]t 


^J2{k—l)n J 


J 2(fc —l)7r 


p2k7T 

< (T^ / cos^{ujit)dt < C 

J 2{k—l)7r 


for some constant C. Since {X^^ are square integrable and {xj^^ are 

bounded in X,F) they are uniformly square integrable. 

As any Levy process is by the Levy-Ito decomposition the sum of a finite variation Levy 
process and an independent square integrable martingale with moments of all orders, we 
obtain the claimed uniform squre integrability for all driving Levy processes. 

Likewise one shows that converges in distribution to J\f ^0, for all 6 G So 

the Cramer-Wold device concludes. 
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Therefore 


(cu^)Yr~"' (wi), zr' (co^)Y^' (w*) 




'•( 2 *) ! 


.( 21 ), 


1 T 


of Theorem 


3.7 


- i=l,...,d 


converges in distribution to 


M ^0, ^I 2 dx 2 d^ and thus using Lemma 3.4 and equations (211, (22) [Y'T{ujj)]J^^ ^ con¬ 
verges to AA (0, ), where B is dehned above. Repeating the reasoning from the proof 


we obtain that ^ converges to a vector of independent, 

for j = □ 


exponentially distributed random variables with Exp 


b{iujj) 


Note that Theorem |3.6| is basically a special case of Theorem 3.8 . However, the proof in 
the case of several frequencies is much more complicated and a more elementary reasoning 
was also presented. 

The limiting result is the analogue of the one for discrete time ARMA models. (See e.g. 
Chapter 10.) 


3.2. Numerical Approximation of Integrals and Limiting Behaviour of the Trun¬ 
cated Pathwise Fourier Transform Based on Non-equidistant Discrete Grids. 

In this section we deal with the numerical approximation of the integral 


(23) Tt{Y){u)^ Y{t)e 


—iujt 


dt. 


Our aim is to describe conditions under which we are able to calculate numerically the 
truncated Fourier transform of a CARMA process based on non-equidistant observations. 
The main result is the following: 

Theorem 3.9. Let X and Y be processes g ive n by the state-space representation ([^ and 
(|^. Suppose that Assumptions 2.1 2.2, 2.t\ and 2.4 are satisfied. Assume that F: M —)• 

be a twice continuously differentiable function with ||F^^||oo < oo. Let 

V / i=0,...,N{T)-2 

(T) (T) 

be a partition of the interval [a, b] with Xq = a and = b and let h^a_ffT) = 

maXj=o,...,7V(r)-i “4"^^ •• 


(24) 


a, 


{N{T)) 

0 


ST) 

'I 


ST) 




a 


(iV(T)) 

N{T)-1 


(T) 


(T) 

‘^N{T)-2 


F (X 


(T) 

"v(r)-i J ’ 


(T) _ (T) 

(25) of = fi±Lj^F (xf)) , j = 1,..., N{T) - 2. 
Then there exist positive constants Ci , C 2 such that 


E 


^ _ f Yit)F{t)dt 

j=o 


2-1 


<CffC2 + T)N{Tfhi,ffT) 


and thus iflimT^ooTN{T)‘^h^g^^{T) = 0, then 


lim 

T^oo 


^( 0-1 ^ ^ .b 

^ (xP^ - j Y{t)F{t)dt 

j=o 


= 0 . 


L 2 
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We begin by establishing an error bound of the trapezoidal method for non-equidistant 
data. For a very accessible approach of quadrature rules procedures we refer to |22| . 
Recall the basic properties of the trapezoidal rule: 

Lemma 3.10. Let /: [a, 6] —)• M 6e a twice continuously differentiable function. Write 

(26) J f{x)dx =’^-^[f{a) + f{b)]+E'^if). 

Then 

(27) \E^{f)\ < sup |/"(x)|. 

For the composite trapezoidal rule for an equidistant grid a < a + {b — a) ^ < ...,a+{b — 
a}^ <■■■< b we have 

(28) ^ mdx = 

Then 

(29) \El{f)\ < f sup |/"(x)|. 

12 n"‘ xela,b] 

A proof can be found e.g. in |22| . 

Now we are going to formulate a version of the trapezoidal rule for non-equidistant points. 
We assume that we have some control on the maximal distance between observations. 

Lemma 3.11. Let a = xo < xi < ■ ■ ■ < xjv_i < xn = b be an arbitrary partition of the 
interval [a, b] and assume that f: [a, 5] —)> M zs a twice continuously differentiable function. 
Put /imax = maxj=o,...,iV- 2 (a;i+i - Xj). Then 

Af-l 

f{x)dx = ^^ [f{xj) + f{xj+i)] + E^{f), 

i=0 

where \E'^{f)\ < iV||/"||oo%^- 



n-l 


2 = 1 


/(«) + 2^/(a-F(6-a)-j f{b) 


+ <(/)• 


Proof 
Let us write 


Af-l 


[a, b]= (j [xj,xj+i], Ij := [xj,Xj+i] 
j=0 


and apply Lemma 3.10 for each interval Lj. Therefore 

Xj+l — Xj 


fXj+i 

/ f{x)dx = 
J X-i 


[f{xj) + f{xj+i)] + E. (/), 


with 


12 

For each z = 0,1,..., — 1 we have 


i?f(/)| < sup \f"{x)\. 

X&[xj,Xj + l\ 


sup \f"{x)\ < sup |/"(x)| =: ||/"| 

xG[xj ,Xj^i] xG[a,b] 


oo • 
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Therefore 


N-l 


f{x)dx = ^ 


Xj + l Xj 


i=0 


where 


\E^{f)\ = 


N-l 




1=0 


N-l 


< iiriioo E 


[f{xj) + f{xj+i)]+E^{f), 


{Xj + l - Xjf 


i=0 


12 


N-l 




i=0 


12 


This completes the proof. □ 

We use some results and ideas from |10] . The aim is to hud an approximation similar to 
Proposition 5.4 of (lO) of the integral appearing in the truncated Fourier transform in the 
case that the observations of the process Y are given on a non-equidistant grid. Let 


N-l 


^ [o,r] 


/ = E 

j=0 


+ fixj+i)] 


be the trapezoidal rule discussed in Lemma 3.11 
stochastic integrals from llOl. 


Recall hrst the Fubini type theorem for 


Lemma 3.12. | 10 [ Theorem 2.4] Let [a,b] C M. be a bounded interval and {L{t))t>o be a 
Levy process with finite second moments. Assume that F: [a, b]xM. ^ is a bounded func¬ 
tion 13{[a, 5])(8)i3([—s, f\)-measurable for all s,t £ (0, oo) and the family {u i—)• F{s, u)}u^[a,b] 
is uniformly absolutely integrable and uniformly converges to zero as |n| —)• 0. Then 


(30) 



F{s, u)dL{u)ds 


F{s, u)dsdL{u) 


a.s. 


In the paper |10] the assumption about measurability in the statement of the theorem is 
not explicitely stated. However, an inspection of their proof combined with results from 
|23| shows that the precise statement has to be in the above form. 


Secondly, note that for non-equidistant data the corresponding error estimation 10, Propo¬ 
sition A. 6 ] has the following form: 


Proposition 3.13. Let [a, 6 ] C M 6 e a compact interval and use the notation of Lemma 

Em 


(1) ///: [a, 6] —)■ M zs twice continuously differentiable, then 

J a 

(2) If F: [a, b] —?• is twice continuously differentiable, then 

J a 

where || • || denotes the Euclidean norm in M'^. 

Here ||F"||oo := maxi=p..,,rfsupi.gja^b] ||F"(ti)||. 


< VdN\\F 




//11 ' ^max 

oo ■ 


12 


<A^II/"lloo^. 
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Put 


g{s)f{s)ds. 


Proof of Theorem 3.9 Assume that we have observed the process Y on the grid 
('T) . ... 

1 ^ -A 

Ar(r)-1 


0 = < ■ ■ ■ < X 1 = T. We have 


( 31 ) T^It]FY= 


j=0 


where {j = 0,..., N{T) — 1) are the coefficients given by (24) and (25). 

Observe that /(a) = f f(s)da(s)ds. Moreover, for all j = 0,... ,N{T) — 1 we know that 
xf"'^ G [0, T], therefore for all u G [0, T] and for all j = 0,..., N{T) — 1 we have 




0,x 


(T) 


(u). 


Thus by (31) we have 

Ar(r)-1 


E «r 

j=0 

N(T)-1 

a 


N{T)-1 


^T) 




mT)H,(jp\ j2 af r b^AYYedL'iu) 

j=0 

(T) 




3=0 


(u) + f 

Jo 


^ - X A.{xYYu 


h^e Yi -VedL*{u) 


aO N(T)-1 

= E ^^)edL*{n) 

• 1 -“ j=0 

f (T) \ 

j=0 

rT^Y)-^ 3 (T) \ . . 

Jo j=o 

rO 

af^'^'^hpT)is)h'^e^^^-YedsdL*{u) 

3=0 

rp p'Y' ^ iY') — 


' —oo Jo 


+ 


' 0 J u 


Y^ ^^edsdL*{u) 


3=0 

rp rp N (T') — 1 


Y^ ^^edsdL*{u). 

4 —oo 4max{0,u} j_Q ^ 

Thus using the representation ^ and the Fubini-type Theorem 3.12 we have 

pT pT ps 

/ F{s)Y{s)ds= / F{s) / b^e^(*-“)edL*(u)ds 

Jo Jo J —oo 
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fT fT 

= / F{s)h'^e^^^-^'^edsdL*iu) 

J— oo t/max{0,w} 
rO pT 

= / F{s)h'^e^^^-^^edsdL*{u) 

J —OO J 0 

pT pT 

+ / F(s)b'^e^("-“)edsdL*(n). 

Jo J u 


Thus 


_ rpN rp-y _ 

^FY ~ ^ 


F{s)Y{s)ds = af^'^^h^,p{s)-F{s)Y{^)ds 

^ -T(s) j b^e^(®"“)edsdL*(u) 


,0 rT (N{T)-1 


l-ocJO \ 

rp rp / N(T'j — 1 


+ J J i ^ — T(s) I edsdL*(u). 


Let us denote 


j=0 


T (N(T)-1 


(32) r*'^^(u) : = J I ^ — F(s) I b^e^^* “^eds, u < 0 , 


j=0 
/N(T)-1 

(33) G^^\u)-.= ^ 

\ i=o 


1n(t)). 


aj ' 5^(t){s) — F{s) I b^e^^* '^^eds, uG [0,T]. 


By Assumption |2.3| we know that there exist positive constants a, f3 such that 
(34) II exp(At)|| </3exp(—at). 

Note that by Lemma [3.11| and Proposition 3.13 we have 

/V(T)-1 \ 

( a^^^'^^h^(T)is) - F{s)\h'^e^^''~^^eds 


fUQ 


N{T)-1 ^ 

(T)(s)b^e^(^^-“)e- / F(s)b^e^("-“)eds 


E 

j=0 


Oi) 


luo 


< Vdiv(r)||T"(u)||oo%^ 

with F{s) = F(s)b'^e'^^^“'“^e and no G [0,r]. If uq = 0, then there exist 5 > 0 and D > 0 
such that ||T"(n)||oo < Zlexp(an) for n < 0. Therefore there exists a constant Di > 0 
such that 

||r('^)(n)b. < VdiV(r)||T|"_^]||oo%^ < ZliiV(T)hLx(r)exp(5n), n < 0. 

If now no = n is any element of [0,T], then there exists D > 0 such that ||T"|[„'j^] ||oo < D 
for n G [0,T]. Therefore there exists a constant D 2 > 0 such that for n G [0,T] we have 

||GW(n)b. < VdiV(r)||F|f„^^]||oo%^ < D2N{T)hl,^{T). 
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By the Ito isometry 

r-0 


T^^\u)dL*{u) 


= E 


L 2 


= E 


= a 



f T^^\u)dL*{u'^ 

’ / /-o \ ^ / 

rO \ 

^T^^\u)dL*{u)j 

1 T^^\u)dL*{u)j 

( 1 ° [r(^\u)fr(^Hu)du^ 


J —OO 



fO 

<(7^ [DiN{T)h^gj^{T) exp{au)) du 
J —OO 


' —OO 

,2 Ar^T-N2i,6 


= a^N(Tyh^^^(T)D( / exp{2au)du 


= DrNiTfhl^yT), 

where > 0 is a constant. In a similar way we obtain 

rT 


G^^\u)dL{u) 


= E 


L 2 


= E 


^ G^^\u)dL{u)^ G^^\u)dL{u 

^ G^^yu)dL{u)^ G^^\u)dL{u] 


= a^^J [G^^Hu)fG^^\u)du 

<a^ r \\G^^\u)\\l,du<a^ f {D2N{T)hl,yT)f du 
Jo Jo 

= DGN{TfThl,yT) 


for some constant Dq > 0. Therefore 

rO 


\E 


,T,Nu2 


FY IIl2 ^ 


< 2 


T^^\u)dL*{u] 


+ 


L 2 


rT 2 

G^^\u)dL{u 


L 2 


thus 


Il 4 y Ili 2 < G, {{G2+T)N{Tfhl,yT)f, 

where Gi, G 2 are positive constants. If limy^oo 'EN{Tyh^^{T) = 0 , then limr_^oo WEpy IIl 2 
0. This completes the proof. □ 


Now we are going to apply Theorem 3.9 to hnd a numerical approximation of the trun¬ 


cated Fourier transform. Using the notation of Theorem 3.9 we denote the trapezoidal 
approximation of 

Et{Y){u;) = ^ £ Y{t)e-^‘^^dt 


N-l 


Tt(Y)(u) = 4 e “5"’!- (4’). 
^ j=o 


by Tt{Y){u), i.e. 
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where the grid points 


^ y j=0,...,N{T)-l 


are given as in Theorem 


3.9 


and 


(T) (T) 

(N{T)) _ -Xq / (t)\ (iV(T)) _ ®JV(T)-1 

“o ~ 2 ~ 


(T) 

''N{T)-2 


F (X 


(T) 

'^Af(T)-l I ’ 


a 


JT) _JT) 
{N{T)) _ ^j+1 ^j-1 


F (xP 


j = l,...,NiT)-2. 


with F{x) = e 

Theorem 3.14. Let X and Y be processes given by the state-space representation ([^ and 
PI). Suppose that Assumptions 2.1, 2.2 2.3 and 2.4 are satisfied and that the process Y T 

'—' -' -' - - (T) (T) 

observed at not necessarily equidistant points Q = Xq < x\ < ■ ■ ■ < x^ 


hrai,y.{T) := maXj=o,...,Ar(T)-2 pP^ - xpy If 


N(T)-1 


IS 

= T. Let 


then 


and thus also 


hm iV(r)hL.(T) = 0 , 

T^oo 


hm \\TT{Y){u)-FT{Y){u)h 2=0 
T^oo 


lim [Tt{Y){u!) - Ft{Y){uj)] = 0. 
T—>-oo 


Proof We identify C with in the canonical way. Applying Theorem 3.9 for 4 = 2, 
F{t) = [cos(a;t), — sin(wt)]'^ we get 

2-1 


E 


V(T)-1 
^ ^(NiT)) 

j=0 


Y (x 


SF 


Y{t) 


cos(wt) 
— sin(t<;t) 


dt 


Dividing both sides by T > 0 we obtain 


E 


N(T)-1 (N(T)) 

j 


E 

j=0 


Vt 


Y X 




1 

-vfl 




2- 

o 

dt 


— sinfwt) 



<Ci(C2 + r)iv(T)2hL.(T). 


< ^+Ciiv(r)2hL.(T). 


Passing to the limit with T —)• oo and using the assumption limj’-^.oo N{T)h^g^^{T) = 0 we 
get the assertion. □ 

Now we are going to state the central limit theorem for the truncated Fourier transform: 

Theorem 3.15. Let X and Y be processes given by the state-space representation (j^ and 

p|). Suppose that Assumptions 2.1 2.2 2.3 and 2.4 are satisfied and the process Y P 

^ - - - iTj (T) 

observed at not necessarily equidistant points 0 = x^ ' < x\ < 


< X 


hma.x{T) := niaXj=o,...,A''(T)- 2 ( 3 ;j+i Let ap^ be defined as in Theorem 

that 


F) — j' 

Assume 


3.9 


hm N{T)hi,,{T) = 0 . 
T^oo 


PutY = ^ 


b{iuj) 

a{iio) 


I 2 x 2 - then 


d — lim 
T^oo 


■3?(rrT(w)) 

^{TtY{oj)) 


= W(0,S), 
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d- lim (3f?(7^y(w))2+ 9(7^y(a;))2) = Exp 


T^oo 

If u = 0, then 

and 


a 


b{iu:) 


a{ioj) 


d - lim TtY{0) = AA f 0, f ® ] ^2 

T —^oo \ 


d — lim 
T-s-oo 


a(0)7^y(0) 


6 ( 0 ) 


/(I). 


Clearly, an analogous statement using Theorem |3.8| holds for the joint distribution when 
the truncated Fourier transform is taken at different frequencies. 

Proof of Theorem 13.151 Put 

Z{T) := r e-^‘^^dL{t) 

and consider the following two-dimensional random vectors: 


Zn :— 


MziT)) 
Q(Z(T)) ’ 


:= 






5?(rry(w)) 

Q(TTr(uj)) 

Observe that it is enough to consider the above limits for T = n. By Lemma [3.4| we know 
that 


From Theorem 3.6 we get 
Therefore 


By Theorem 3.14 we have 
Therefore, 

In the same way we obtain 


- lim llUn - Znll = 0. 


d— lim Z„=AA(0, S). 


d— lim Un=AA(0, S). 
n^oo 


P - lim (V„ - Un) = 0. 

n^oo 

d- lim =AA(0,S). 


d — lim |Z|2 = Exp cj^ 


b{iu) 


a{iu)) 


In order to obtain the assertion for w = 0 we repeat the above resonings applying Theorem 
3.51 instead of Theorem 13.41 □ 


4. Illustrative simulations 

We now turn to a numerical illustration of the theoretical convergence results given in 


Section 3.2 We are looking at simulations of CARMA processes and their numerically 
approximated truncated Fourier transform over different time horizons and maximal grid 
widths. To illustrate the convergence to the asymptotic normal distribution we shall look 
at several frequencies and different driving Levy processes, standard Brownian motion, 
a Variance Gamma process and a “two sided Poisson process”. Of course, the truncated 
Fourier transform of (Tj) is obtained using the trapezoidal rule based on non-equidistant ob¬ 
servations of the CARMA process (Tj) given on the interval [0, T]. On an interval [0, T] we 
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generate a non-equidistant grid in the following way: we fix the maximal distance hma.x{T) 
between elements of the grid and from each interval \i ■ (i + 1) • ^hinax(^)) for 

i = 0,1,..., — 1 we draw a nnmber according to the nniform distribntion. This resnlts 

in a non-eqnidistant grid with the nnmber of points being N(T) = 2T//imax + 1- 

For onr simnlations we nsed the R Project for Statistical Compnting. For the simnlation 
we first generate the non-eqnidistant grid by the above procednre and then join it with a 
regnlar grid of mesh 0.001, which is still on average hve times hner than the non-eqnidistant 
grid of the largest time horizon considered. On this joint grid the CARMA process Y is 
simnlated with a standard Enler scheme for the state space representation. Afterwards 
only the simnlated values at the times of the original non-eqnidistant grid are nsed to 
compnte the approximation of the trnncated Fonrier transform with the trapezoidal rnle. 
In all cases we simnlate 2000 independent paths of the CARMA process and compnte the 
associated valnes of the trnncated Fonrier transform at the following freqnencies: 


[wi, W2,W3,a;4] = [0, 0.1, 1, 10]. 


For the non-zero freqnencies real and imaginary part have to be considered separately. 
However, in the following we look only at the real parts as the behaviour of the imaginary 
parts is most similar. Mainly, the resnlts are presented via QQ-plots where the theoretical 


valnes follow the (limiting) law described in Theorem 3.15 


We are going to consider CARMA processes with the following antoregressive and moving 
average orders: (p, q) = (1, 0), i.e. an Ornstein-Uhlenbeck type process, and (p, q) = (2,1). 
For the time horizon T and the maximnm distance of the non-eqnidistant observation times 
we consider the pairs (T = 10, /imax = 0.1), (T = 50, /imax = 0.05) and (T = 100, /imax = 
0 . 01 ). 


For each case we consider three different driving Levy noises: standard Brownian Motion, 
a Variance Gamma process and a “two sided Poisson process”. For the dehnition and 
properties of the Variance Gamma process we refer to (l^ and references therein. We 
constrnct the process in the following way: V) = G^ — G^, where Gj and G^ are independent 
Gamma processes with shape parameter 1 and scale parameter 4. Likewise the “two sided 
Poisson process” is the difference of two independent Poisson processes with rate 10, i.e. a 
componnd Poisson process with rate 20 and jnmps -|-1 and —1 both with probability 1/2. 


Example 4.1. We consider the CAR(l) model. Then A = —ai and 


a{z) := z + ai, b(z) = bo. 


So the spectral density is 


/(w) 


cr^ b{iuj) 
In a{iu}) 


bl 

27r -|- af 


For the simnlations we take [bo, ai] = [1, 2]. 

QQ-plots showing the resnlts for 2000 simnlated paths for the fonr different freqnencies 
and three different combinations of time horizon and maximnm grid width can be found 
in Fignres Eli and for the driving Levy process being a standard Brownian motion, a 
Variance Gamma and a two-sided Poisson process, respectively. 


Example 4.2. We consider the GARMA(2,1) model. We have 


0 1 

— 02 —Ol 



e := 


0 

1 ’ 




A := 
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The autoregressive and moving-average polynomials are of the form 

a{z) = z'^ + aiz + 02 , b{z) = z + bo. 


We have 
b{iu) 
a{iLo) 


ioj + bp f( ) - ^ 

{iioy + {iio)ai + a2' 27r 


b{iuj) 


a{ioj) 


a 


bl + 


OJ 


2tt ui'^ + (of — 202 ) 0 ;^ -|- a| 


For the simulation procedure we take [ 60 , fei, oi, 02 ] = [1,1,1, 2], 


QQ-plots showing the results for 2000 simulated paths for the four different frequencies 
and three different combinations of time horizon and maximum grid width can be found 
in Figures]^ and for the driving Levy process being a standard Brownian motion, a 
Variance Gamma and a two-sided Poisson process, respectively. Likewise, Figures [^and 
l^show corresponding histograms. 


The simulation results seem to indicate the following. 

In Figure we notice at hrst a pretty good ht of the empirical quantiles from the simula¬ 
tions with the theoretical ones of the asymptotic distribution across all time horizons and 
frequencies. Looking more carefully, the ht in the tails clearly improves when the time 
horizon/hneness of the grid increases, but it is never bad. For the longest time horizon 
and hnest grid the ht is clearly very good. Of course, it should not be forgotten that 
in this case the distribution of the (trapezoidal approximation of the) truncated Fourier 
transform is always exactly Gaussian and not only asymptotically. When looking across 
the non-zero frequencies one notes that for the shortest time horizon the quantiles for the 
smallest frequency 0.1 appear to lie on a line which is somewhat different from the line of 
the theoretical quantiles. This indicates that the quantiles of the simulated paths come 
from a normal distribution, but one with a different variance then the asymptotic one. It 
is no surprise that this occurs for the lowest frequency and the smallest time interval, as 
for low frequencies one observes - regardless of the hneness of the sampling - the fewest 
full cycles over a time interval of hxed length. For this combination of time horizon and 
frequency we see only one full cycle. 

Turning to Figure we hrst notice that the ht in the tails improves again clearly with 
increasing (T, l//imax)- Especially, for the highest (T, l/hmax) one sees that the ht in 
the tails is a bit worse now for a driving Variance Gamma process compared with the 
driving Brownian motion in Figure Of course, now the simulated values are indeed 
only asymptotically following a Gaussian distribution. Looking at the different non-zero 
frequencies one again sees that the ht improves for the higher frequencies. Most notably 
for the lowest frequency one sees for the smallest T = 10 again that the points do seem 
to he on a straight line in the normal QQ-plot, but one with a different slope than for 
the theoretical quantiles. Hence, the variance is clearly different from the asmpytotic one. 
Obviously, this effect is now more pronounced than in the case of the driving Brownian 
motion. 

Moving on to the case of the driving process being a two-sided Poisson process in Figure]^ 
we hrst of all note that again the ht in particular in the tails clearly improves with increasing 
(r, l//ijnax)- For the highest (T, l/hmax) the simulated and theoretical asymptotic quantiles 
agree again extremely well. Again it is certainly a bit worse than in the case of a driving 
Brownian motion, but it seems to be very similar to the Variance Gamma case, although 
maybe for frequency 0 the agreement of the quantiles is slightly worse. Turning to the 
behaviour across non-zero frequencies, we see again that the quantiles are closer for higher 
frequencies and that for the smallest non-zero frequency and time horizon the empirical 
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Normal Q-Q Plot 


Normal Q-Q Plot 


Normal Q-Q Plot 



Theoretical Quantiles 

(a) a; = 0,r=10 

Normal Q-Q Plot 


Theoretical Quantiles 


(d) w = 0.1,r = 10 

Normal Q-Q Plot 


(g) w = 1,T = 10 

Normal Q-Q Plot 


(j) w = io,r = 10 






Theoretical Quantiles 

(b) u! = 0,T = 50 

Normal Q-Q Plot 



Theoretical Quantiles 


(e) w = 0.1, T = 50 

Normal Q-Q Plot 



(h) w= 1,T = 50 

Normal Q-Q Plot 



(k) w = 10,T = 50 



(c) a; = 0,r= 100 

Normal Q-Q Plot 



Theoretical Quantiles 


(f) w = 0.1, r = 100 

Normal Q-Q Plot 



(i) w = i,r = 100 

Normal Q-Q Plot 



(l) w = 10,T= 100 


Figure 1. Normal QQ plots for the real part of the truncated Fourier 
transform of the Ornstein-Uhlenbeck type process driven by stan¬ 
dard Brownian Motion for the frequencies 0,0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns). The theoretical quantiles are coming from the (limiting) law 


described in Theorem 3.15 
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Normal Q-Q Plot 


Normal Q-Q Plot 


Normal Q-Q Plot 



Theoretical Quantiles 

(a) w = 0,T = 10 

Normal Q-Q Plot 


(d) oj = 0.1,T= 10 

Normal Q-Q Plot 


Theoretical Quantiles 

(g) uj = 1,T= 10 

Normal Q-Q Plot 


(j) w = io,r = 10 






Theoretical Quantiles 

(b) w = 0,r = 50 

Normal Q-Q Plot 



(e) a; = 0.1,T = 50 

Normal Q-Q Plot 



Theoretical Quantiles 

(h) w = 1,T = 50 

Normal Q-Q Plot 



(k) w = 10,T = 50 



(c) w = 0,T = 100 

Normal Q-Q Plot 



(f) w = 0.1,T=100 

Normal Q-Q Plot 



Theoretical Quantiles 

(l) w = 1,T= 100 

Normal Q-Q Plot 



(l) w = 10,T = 100 


Figure 2. Normal QQ plots for the real part of the truncated Fourier 
transform of the Ornstein-Uhlenbeck type process driven by a Vari¬ 
ance Gamma process for the frequencies 0,0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns). The theoretical quantiles are coming from the (limiting) law 


described in Theorem 3.15 
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Normal Q-Q Plot 


Normal Q-Q Plot 


Normal Q-Q Plot 



Theoretical Quantiles 

(a) a; = 0,r=10 

Normal Q-Q Plot 


(d) a; = 0.1,r= 10 

Normal Q-Q Plot 


Theoretical Quantiles 

(g) w = 1,T = 10 

Normal Q-Q Plot 


(j) w = io,r = 10 






Theoretical Quantiles 

(b) w = 0,r = 50 

Normal Q-Q Plot 



(e) a; = 0.1,T = 50 

Normal Q-Q Plot 



Theoretical Quantiles 

(h) 1,T = 50 

Normal Q-Q Plot 



(k) w = 10,T = 50 



(c) a; = 0,r= 100 

Normal Q-Q Plot 



(f) w = 0.1, r = 100 

Normal Q-Q Plot 



Theoretical Quantiles 

(i) w = i,r = 100 

Normal Q-Q Plot 



(l) w = 10,T= 100 


Figure 3. Normal QQ plots for the real part of the truncated Fourier 
transform of the Ornstein-Uhlenbeck type process driven by a two- 
sided Poisson process for the frequencies 0,0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns). The theoretical quantiles are coming from the (limiting) law 


described in Theorem 3.15 
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Normal Q-Q Plot 


Normal Q-Q Plot 


Normal Q-Q Plot 



Theoretical Quantiles 

(a) w = 0,T = 10 

Normal Q-Q Plot 


(d) oj = 0.1,T= 10 

Normal Q-Q Plot 


Theoretical Quantiles 

(g) uj = 1,T= 10 

Normal Q-Q Plot 


(j) w = io,r = 10 






Theoretical Quantiles 

(b) w = 0,r = 50 

Normal Q-Q Plot 



(e) a; = 0.1,T = 50 

Normal Q-Q Plot 



Theoretical Quantiles 

(h) w = 1,T = 50 

Normal Q-Q Plot 



(k) w = 10,T = 50 



(c) w = 0,T = 100 

Normal Q-Q Plot 



(f) w = 0.1,T=100 

Normal Q-Q Plot 



Theoretical Quantiles 

(l) w = 1,T= 100 

Normal Q-Q Plot 



(l) w = 10,T = 100 


Figure 4. Normal QQ plots for the real part of the truncated Fourier 
transform of the simulated CARMA(2,1) processes driven by stan¬ 
dard Brownian Motion for the frequencies 0,0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns). The theoretical quantiles are coming from the (limiting) law 


described in Theorem 3.15 
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Histogram of data Histogram of data Histogram of data 



(A)w = o,r = io (b)w = o,t = 50 (c)w = o,r = ioo 

Histogram of data Histogram of data Histogram of data 



(d) w = 0.1,T= 10 (e) u; = 0.1,T = 50 (f) w = 0.1, T = 100 

Histogram of data Histogram of data Histogram of data 



(G)a; = l,T = 10 (h)w=1,T = 50 (l) w = 1, T = 100 

Histogram of data Histogram of data Histogram of data 



(j)u; = 10,T = 10 (k)w = 10,T = 50 (l) w = 10, T = 100 

Figure 5. Histograms and limiting density for the real part of the trun¬ 
cated Fourier transform of the simulated C ARM A (2,1) processes driven by 
standard Brownian Motion for the frequencies 0, 0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns) 
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Normal Q-Q Plot 


Normal Q-Q Plot 


Normal Q-Q Plot 



Theoretical Quantiles 

(a) w = 0,T = 10 

Normal Q-Q Plot 


(d) oj = 0.1,T= 10 

Normal Q-Q Plot 


Theoretical Quantiles 

(g) uj = 1,T= 10 

Normal Q-Q Plot 


(j) w = io,r = 10 






Theoretical Quantiles 

(b) w = 0,r = 50 

Normal Q-Q Plot 



(e) a; = 0.1,T = 50 

Normal Q-Q Plot 



Theoretical Quantiles 

(h) w = 1,T = 50 

Normal Q-Q Plot 



(k) w = 10,T = 50 



(c) w = 0,T = 100 

Normal Q-Q Plot 



(f) w = 0.1,T=100 

Normal Q-Q Plot 



Theoretical Quantiles 

(l) w = 1,T= 100 

Normal Q-Q Plot 



(l) w = 10,T = 100 


Figure 6. Normal QQ plots for the real part of the truncated Fourier 
transform of the simulated CARMA(2,1) processes driven by a Vari¬ 
ance Gamma process for the frequencies 0,0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns). The theoretical quantiles are coming from the (limiting) law 


described in Theorem 3.15 
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Histogram of data Histogram of data Histogram of data 



(A)w = o,r = io (b)w = o,t = 50 (c)w = o,r = ioo 

Histogram of data Histogram of data Histogram of data 



(d) w = 0.1,T= 10 (e) u; = 0.1,T = 50 (f) w = 0.1, T = 100 

Histogram of data Histogram of data Histogram of data 



(G)a; = l,T = 10 (h)w=1,T = 50 (l) w = 1, T = 100 

Histogram of data Histogram of data Histogram of data 



(j)u; = 10,T = 10 (k)w = 10,T = 50 (l) w = 10, T = 100 

Figure 7. Histograms and limiting density for the real part of the trun¬ 
cated Fourier transform of the simulated C ARM A (2,1) processes driven by 
a Variance Gamma process for the frequencies 0, 0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns) 
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Normal Q-Q Plot 


Normal Q-Q Plot 


Normal Q-Q Plot 



Theoretical Quantiles 

(a) w = 0,T = 10 

Normal Q-Q Plot 


(d) oj = 0.1,T= 10 

Normal Q-Q Plot 


Theoretical Quantiles 

(g) uj = 1,T= 10 

Normal Q-Q Plot 


(j) w = io,r = 10 






Theoretical Quantiles 

(b) w = 0,r = 50 

Normal Q-Q Plot 



(e) a; = 0.1,T = 50 

Normal Q-Q Plot 



Theoretical Quantiles 

(h) w = 1,T = 50 

Normal Q-Q Plot 



(k) w = 10,T = 50 



(c) w = 0,T = 100 

Normal Q-Q Plot 



(f) w = 0.1,T=100 

Normal Q-Q Plot 



Theoretical Quantiles 

(l) w = 1,T= 100 

Normal Q-Q Plot 



(l) w = 10,T = 100 


Figure 8. Normal QQ plots for the real part of the truncated Fourier 
transform of the simulated C ARM A (2,1) processes driven by a two 
sided Poisson process for the frequencies 0,0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns). The theoretical quantiles are coming from the (limiting) law 


described in Theorem 3.15 
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Histogram of data Histogram of data Histogram of data 



(A)w = o,r = io (b)w = o,t = 50 (c)w = o,r = ioo 

Histogram of data Histogram of data Histogram of data 



(d) w = 0.1,T= 10 (e) u; = 0.1,T = 50 (f) w = 0.1, T = 100 

Histogram of data Histogram of data Histogram of data 



-3 -2 -1 0 1 2 3 4 -4 -3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3 

data data data 


(G)a; = l,T = 10 (h)w=1,T = 50 (l) w = 1, T = 100 


Histogram of data Histogram of data Histogram of data 



(j)u; = 10,T = 10 (k)w = 10,T = 50 (l) w = 10, T = 100 

Figure 9. Histograms and limiting density for the real part of the trun¬ 
cated Fourier transform of the simulated C ARM A (2,1) processes driven by 
a two-sided Poisson process for the frequencies 0,0.1,1,10 (rows) and time 
horizons/maximum non-equidistant grid sizes 10/0.1,50/0.05,100/0.01 
(columns) 
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quantiles seem to be in line with a normal distribution with a somewhat different variance 
compared to the asymptotic one. The size of this effect seems to be rather similar to the 
Variance Gamma case. Interestingly, also at frequency 0 the QQ-plot seems to indicate for 
T = 10 that the empirical quantiles are close to the ones of a normal distribution with a 
slightly different variance than the asymptotic one. 

To summarize the simulation study in the CAR(l)/OU-type case we can clearly conclude 
that the asymptotic distribution result approximates the finite-sample distribution of the 
trapezoidal approximation of the truncated Fourier transform in our simulations very well 
and that the convergence to the asymptotic distribution is fast. For small frequencies, 
especially when one has about one full cycle or less over the time horizon considered, one has 
to be careful, as then the distribution tends to be somewhat different from the asymptotic 
one for good reasons. This effect seems to be more pronounced when one considers a 
Levy process with jumps compared to a Brownian motion. In general the quality of the 
approximation of the simulated quantiles by the asymptotic ones is somewhat better in 
the case of a Brownian motion than in a pure jump process. Comparing the driving jump 
processes, the finite activity rather discrete two-sided Poisson process with the infinite 
activity Variance Gamma process, we do not see any significant differences. It should 
be noted that both jump processes are, however, light-tailed in the sense that they have 
exponential moments. It would not be surprising if this picture changes when considering a 
really heavily tailed driving Levy process. Note that our theoretical results are valid also in 
rather heavily-tailed cases. For the asymptotic normality of the (trapezoidal approximation 
of the) truncated Fourier transform we only needed finite second moments. 

Turning to the simulations of CARMA(2,1) processes, most of the findings of the CAR(l)/ 
OU case remain valid, so we only point out the differences. In the case of a driving Brownian 
motion, depicted in Figure]^ the only difference seems to be that for T = 10 and u = 0.1 
the empirical quantiles are now appearing to lie on a line farther away from the theoretical 
quantiles which implies that in the G ARM A (2,1) case the variance in the simulations is 
clearly farther away from the asymptotic one than in the OU case. The same applies for 
the Variance Gamma case of Figure and the two-sided Poisson case of Figure On 
top of the QQ plots we now also provide histograms in Figures i0 and 1^ respectively, 
together with plots of the limiting normal density. To us it seems very hard to see the 
convergence to normality with increasing T in the histograms, which reflects the fact that 
it is essentially the tails which need to converge and they are much clearer visible in the 
QQ plots than in histograms. It is also not easy to see in them that for cj = 0.1, T = 10 
the variance of the simulated values is different from the asymptotic theoretical one. The 
only thing one notices is that for w = 0.1, T = 10 the histogram routine of R tends to 
use very different bins than in all the other cases. Note that all histograms were obtained 
using the default parameters of the hist function in R, so the binning was done by the 
standard automatic selection to give “nice” histograms. Hence, from our simulations of 
GARMA(2,1) processes we can conclude that the orders of the GARMA processes and the 
particular autoregressive and moving average parameters appear not to really matter for 
the (qualitative) behaviour of the (trapezoidal approximation of the) truncated Fourier 
transform. 


5. Conclusion and Outlook 

We have obtained an asymptotic normality result for the (trapezoidal approximation of 
the) truncated Fourier transform under essentially minimal assumptions (i.e. second mo¬ 
ments) and seen via a simulation study that this result approximates the finite sample 
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behaviour very well, unless the frequency is too low compared to the length of the con¬ 
sidered time interval. This suggests clearly that it should be very promising to develop 
statistical inference techniques for non-equidistantly sampled CARMA processes by con¬ 
sidering continuous observation techniques and using numerical approximation schemes to 
compute the quantities of interest based on the observed non-equidistant data. The ap¬ 
propriate set-up to get convergence and asymptotic distribution results is to send the time 
horizon to infinity and to send at the same time the maximum distance of observation time 
points to zero. 

Based on our results in this paper it seems natural to locally smooth the trapezoidal ap¬ 
proximation of the truncated Fourier transform to get consistent estimators of the spectral 
density and to use it in a Whittle type estimator for the AR and MA parameters. Consid¬ 
ering this is beyond the scope of the present paper. 
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